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ABSTRACT 

Gamma Ray Bursts (GRBs) are the most energetic sources in the universe 
and among the farthest known astrophysical sources. These features make them 
appeahng candidates as standard candles for cosmological apphcations so that 
studying the physical mechanisms for the origin of the emission and correlations 
among their observable properties is an interesting task. We consider here the 
luminosity L\ - break time T* (hereafter LT) correlation and investigate whether 
there are systematics induced by selection effects or redshift dependent calibra- 
tion. We perform this analysis both for the full sample of 77 GRBs with known 
redshift and for the subsample of GRBs having canonical X-ray light curves, 
hereafter called f/0095 sample. We do not find any systematic bias thus con- 
firming the existence of physical GRB subclasses revealed by tight correlations 
of their afterglow properties. Furthermore, we study the possibility of applying 
the LT correlation as a redshift estimator both for the full distribution and for 
the canonical lightcurves. The large uncertainties and the non negligible intrin- 
sic scatter make the results not so encouraging, but there are nevertheless some 
hints motivating a further analysis with an increased f/0095 sample. 
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Introduction 



The high fluence values (from 10 ^ to 10 ^ erg/cm^) and the huge isotropic energy emit- 
ted 10^° - lO^^erg) at the peak in a remarkably short prompt emission phase make GRBs 
the most violent and energetic astrophysical phenomena. Fifty years after their discovery 
in the '60s by the Vela satellites, the nature of GRBs is still unclear. Notwithstanding the 
variety of their peculiarities, some common features may be identified by looking at their 
light curves. GRBs have been traditionally classified as short (Tgo < 2s) and long (Tgo > 2s), 
although some recent studies (see, e.g., Norris & Bonnell 2006) have revealed the existence 
of an intermediate class (IC) thus asking for a revision of this criterium. Consequently, long 
GRBs have now been divided into two classes, normal and low luminosity, the latter one 
likely being associated with Supernovae. Here, we concentrate our attention on the class of 
normal long bursts, observed in X-ray with the aim of better clarifying their origin in the 
view of possible systematics. A valid tool in classifying GRBs is pro vided by th e analysis 
of their light curves. The data observed by the Beppo - Sax satellite (jPiro 200l[ ) were rea- 
sonably well fitted by a simple phenomenological power- law expression, f{t,v) oc t^^v^^^ 
with (a,/3) — (1.4,0.9). However, a crucial breakthrough in this field has been represented 
by the launch of the Swift satellite in 2004. The Swift instrumental setup, composed by 
the Burst Alert Telescope (BAT, 15 - 150 keV), the X-Ray Telescope (XRT, 0.3 - 10 keV) 
and the Ultra - Violet/ Optical Telescope (UVOT, 170 — 650 nm), allows a rapid follow-up of 
the afterglows in different wavelengths giving better coverage of the GRB light curve than 
the previous missions. Such data revealed the existence of a more complex phenomenology 
with different slopes and break times thus stressing the inadequacy of a single power- law 
function. A significant step forward has been made by the analysis of the X-ray afterglow 
curves of th e full sample of Swift G RBs showing that they may be fitted by a single analytical 
expression ( IWillingale et al. 20071 ) which we referred to in the following as the W07 model. 



It is worth stressing that finding out a universal feature would allow us to recognize 
if GR Bs are standard ca ndles lookin g for correlations amon g their obs ervables. The Ej^.^ - 
Epeak ( Amati et al. 20091), E^-Ep^gk ( Ghirlanda et al. 2006 ). L-Ep^ak ( Schaefer 2003 ) and 
L-V (IRiechart et al. 200 ll ) correlations are some of the attempts pursued in this direc- 
tion. However, the problem of la rge data scatter in the considered luminosity relations 
( iButler et al. 20091 : IYu et al. 20091) and a possible imp act of detector thresholds on cosmo- 
l ogical standard cand les (IShahmoradi fc Nemiroff 2009[ ) have been discussed controversially 
( jCabrera et al. 20071 ). Within this wide framework, we consider here the LT correlation be- 
tween the break time T* = Ta/ (1 + z) and the luminosity at the break time L*x where z 
is the GRB redshift and with asteriks we refer to the rest frame quantities. Both these 
observables refer to the plateau phase of the W07 model. Dainotti et al. (2008) first 
found that these quantities are not independent, but rather follow the log- linear relation. 
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\ogLx = alogT* + 6, with a and b fixed by the fitting proc edure. The LT correlation has 
then b een confirmed (IGhiselhni et aL 2008t lYamazaki 20091 ) and recently updated with 77 
GRBs (jPainotti et al. 2010l ) leading to the discovery of a new subclass of the afterglows with 
smooth observed X-ray light curves, which are prefentially distributed at higher luminosities 
than the full distribution. 

The plan of the paper is as follows. In Section 2, we review the LT correlation explaining 
how the interested quantities are evaluated and the calibration procedure adopted. Selection 
effects are discussed in Section 3, while the problem of a possible evolution with z of the 
calibration parameters is addressed in Section 4. Section 5 investigates the possibility of 
using the LT correlation as a redshift indicator, while a summary of the results is finally 
given in Section 6. 



2. The LT correlation 

The LT correlation relates the time scale T* and the X-ray luminosity L*x at T^, where 
Ta is defined as the end o f the plateau phase. Having had the confirmation of the existence 



of the above correlation (jPainotti et al. 2010l ). we here try to answer the question: Is it 

affected by systematics ? 

As a preliminary remark, let us remember how the quantities of interest are evaluated. 
The source rest frame luminosity in the Swift XRT bandpass, {Emm, Emax) = (0.3, 10) keV, 
is computed as : 



L*x{Emin, Emax, t) = A'kD\{z) Fx{Emin, Emax, t) ■ K (1) 

where Dl{z) is the GRB luminosity distance at redshift z, Fx is the measured X-ray energy 
flux (in erg/cm^/s) and K is the K - correction. Denoting with f{t) the Swift light curve and 
following Bloom et al. (2001), we get : 



iE:D{i:>^E)dE 



KFx{Em^n, Emax, t) = f {t) X -^(1+^) - ' (3) 

f^™- E^(E)dE 

with $(-E) the differential photon spectrum. We model this term as ^[E) oc E^^'^ oc 
where {P>a,la) are the spectral and photon index, respectively. It is worth stressing that the 
fit of the model $(-E) is performed considering only the spectrum of the plateau phase, 
selected using a filter time fixed as Tg ± gt„ \ the Tg values t ogether with their errorbars, ot^ , 
are derived in the fitting procedure (jWillingale et al. 20071 ). As shown also in previous XRT 
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spectral analysis (INousek et al. 20061 ) . this particular choice of the filter time leads to the 
single power - law function as a better fit than the more commonly assumed Band function 
( iBand et al. 1993[ ). According to the W07 model, the functional expression for f{t) is: 



fit) = Uit) + fait) (3) 

where the first term accounts for the prompt (the index "p") 7 -ray emission and the initial 
X-ray decay, while the second one describes the afterglow (the index "a"). Both components 
are modeled with the same functional form : 



fit) 



Fc exp ( etc - ^ 1 exp 



Tr 



-J 1 for t < Tr 



(4) 



exp - 



t, 



for t > Tr 



t 

JcJ ^ \ t 

where c = p 01 a. The transition from the exponential to the power law occurs at the point 
(Tc, Fc) where the two functional sections have the same value and gradient. The parameter 
etc is the temporal power law decay index and the time tc is the the initial rise time scale. 
We refer to Willingale et al. (2007) for further details on the analysis, while we only remind 
here that a usual fitting of the log if lux) vs log itime) data provides estimates and 
uncertainties on the time parameters (logTp, logT^) and the products (logFpTp, logF^Ta). 

For the afterglow part of the light curve, we have computed values L*x (eq. E]) at the 
time Ta, which marks the end of the plateau phase and the beginning of the last power law 
decay phase. We have considered the following approximation which takes into accounts the 
functional form, /„, of the afterglow component only: 



fiTa)^faiTa) = Faexp\^-^j (5) 

where we set ta = Tp because in most cases the afterglow component is fixed at the transition 
time of the prompt emission, Tp. Actually, we are using Eq.(|5]), instead of ([3]) since the 
contribution of the prompt component is typically smaller than 5%, much lower than the 
statistical uncertainty on faiTa)- Neglecting fpiTa) thus allows to reduce the error on FxiTa) 
without introducing any bias. This latter error is then estimated by simply propagating those 
on Pa, ^ogTa and logF^Ta thus implicitly assuming that their covariance is null. Inserting 
Eqs.(l5]) and ([2]) into Eq.([l]), one then obtains: 
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AnDl{z)Fx 

where Fx = Faexp {—Tp/Ta) is the observed flux at the time T^. 

As a final important remark, we note that the presence of the luminosity distance Dl{z) 
in Eq.([6]) constrains us to adopt a cosmological model to compute L*x. We then use a fiat 
ACDM model so that the luminosity distance reads : 

Dl{z) = ^{1 + z) , (7) 



In agreement with the WMAP seven year results (IKomatsu et al. 20101 ) . we set {QM,h) 



(0.272,0.704) with h the Hubble constant Hq in units of 100 km/s/Mpc. 

2.1. Calibration parameters 

Let us suppose that R and Q are two quantities related by a linear relation 



R = aQ + b 



and denote with Cint the intrinsic scatter around this relation. Calibrating such a relation 
means determining the two coefficients {a,b) and the intrinsic scat ter o"j„(. To this aim, 
we will resort to a Bayesian motivated technique (ID' Agostini 20051 ) thus maximizing the 
likelihood function £(a, b, amt) = exp [— L(a, 6, aint)] with : 



Li a, 6, Gi 



int ) 



■ 2 I 2,22 
.^int + ^ft +« ^Q, 



{R, - aQi - bf 



(9) 



where the sum is over the JV objects in the sample. Note that, actually, this maximization 
is performed in the two parameter space (a, aint) since b may be estimated analytically as : 



Ri — aQi 



cr, 



int 



+ 4. + 



int 



-1 



(10) 



so that we will not consider it anymore as a fit parameter. The above formulae easily 
applies to our case setting R = log L*j^(Ta) and Q = logT^. We estimate the uncertainty on 
log L*x{Ta) by propagating the errors on {Ta,Tp, FaTa, /3a). 
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The Bayesian approach used here also allows us to quantify the uncertainties on the fit 
parameters. To this aim, for a given parameter pi, we first compute the marginalized like- 
lihood Ci{pi) by integrating over the other parameter. The median value for the parameter 
Pi is then found by solving : 



Pi,med rPi 

Ci{pi)dpi = - Ci{pi)dpi . (11) 



2 

Pi,min ^ V' 



i.rmn 



The 68% (95%) confidence range {pi,hPi,h) are then found by solving: 



Ci{pi)dpi = — — / Ci{pi)dpi , (12) 



2 

Pi, I ^ Pi.min 



'Pi,h \ — e 



-'I, max 



i.rmn 



Pi,med ^ P 

with e = 0.68 (0.95) for the 68% (95%) range respectively. 



Ci{pi)dpi = I Ci{pi)dpi , (13) 



3. Threshold selection of the fit error parameter 

Dainotti et al. (2010, hereafter DIO) have recently updated the LT correlation using a 
sample of 77 GRBs with known re dshift and Swift X-ray afterglow light curves. DIO have 
defined a fit error parameter u = ^J'ajT~+~a^, measured in the burst rest frame, to analyze 

how accuracy of fitting the canonical lightcurve, (eq. [3] and H]) to the data influences the 
studied correlations. This definition is used to distinguish the canonical shaped light curves 
from the more irregular ones, perturbed by secondary flares and various non uniformities. 
DIO have then defined a fiducial sample selecting only GRBs with u < 4 and excluding the 
IC objects thus selecting 62 out of the original 77 GRBs. To be consistent with DIO, we 
here still consider only the fiducial sample. As a general remark, we would like to stress that 
the study of a whatever correlation among GRBs observables should involve only physically 
homogenous subsamples thus motivating our exclusion of the IC GRBs because of their 
different properties from the long ones that mainly constitute our sample. In other words, 
with homogenous sample we indicate a subsample of GRBs that have lightcurves well defined, 
in the sense that the Willingale model represents with very good accuracy the parameters 
values representing the afterglow and the plateau. As a consequence this subsample tightly 
obeys to the correlation and from this evidence we infer that the properties of the GRBs in 
this subset are the same. For example, the XRFs are included in the subsample, since they 
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Fig. 1. — Spearman correlation coefficients for the L*x-T* (black squares) and (3a- T* (red 
asterisks) as function of the threshold u value. 



obey to the correlations, giving in this way evidence to the theory according to which they 
have the same progenitor mechanism of the normal long GRBs. 

As a ffi st indicator for the existence of a relation, we use the Spearman correlation 
coefficient p (ISpearman 1904| ) providing a non - parametric measure of the statistical signifi- 
cance of the dependence between two quantities. Fig.[T]shows plt = p{L*x,T*) as a function 
of the threshold Uth value, used to exclude from the fiducial sample GRBs with u > Uth- 
As we can note, the smaller uth is, the larger p is, i.e. the more we are confident that a 
statistically meaningful correlation indeed exists. The same figure also shows a similar anal- 
ysis for pi3T = p{(3a,T*) suggesting that also the slope of the GRB spectrum and the 
break time T* are correlated. It is worth noting that the smaller u, the smaller the error 
on {logT*,logL*x) is. Examining the light curves, we find that small errors are obtained 
for the GRBs that follow better the W07 model. We therefore argue that both the LT and 
Pa - T* correlations are statistically meaningful provided the GRBs in the sample belong to 
the class well described by the W07 model. 



In order to better investigate the impact of the u selection, we fit the LT correlation to 



50.0 



log «,|, 



Fig. 2. — Number of GRBs, calibration parameters (a, b) and intrinsic scatter aint as function 
of the threshold Uth value. See text for a detailed discussion. 

GRBs subsamples obtained by selecting only those objects with u < Uth with uth running 
from 0.095 to 4 in steps of 0.01 (with uth = 4 for the fiducial sample and Uth = 0.095 for 
the canonical ones). The upper left panel in Fig. [2] shows that the number of GRBs in the 
sample obviously increases with Uth, but the price to pay is including GRBs with large errors 
on (logT^, \ogL*x). Such large uncertainties may be due to bad sampling or to a less precise 
determination of the parameters fitted within the W07 model. In both cases, the estimated 
values of the fit parameters (logTp, logT^, logF^Ta) are not reliable so that it is a safer option 
to not include large u GRBs in the analysis of correlations. Our chosen value u = 4 represents 
a compromise between the need to assemble a statistically meaningful sample and avoiding 
uncertain couples (log T^, log L^^^) that can unnecessarily increase the intrinsic scatter. 

As a first interesting result, we find that the intrinsic scatter aint is smaller for smaller 
Uth selected samples, with a sharp drop for Uth = 0.4. The high values of plt and the 
decreasing scatter point towards a scenario where the GRBs most deviating from the LT 
correlation are actually the ones with the less precise determination of the fitted parameters 
consistent with our guess that their estimated values of logT^ and logL*x are not reliable. 

It is worth wondering whether selecting on u biases the calibration parameters {a,b). 
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Fig. 3. — Distributions of (a,6, dj^t) and — InL, defined by Eq.Q, for the fiducial sample 
(62 GRBs with m < 4). 



Upper right and lower left panels in Fig.|2] indeed shows a clear increase of b as Uth gets 
smaller, while the slope a of the correlation remains almost constant at the value a ~ —1.06, 
consistent with the results in DIO that the small u GRBs (referred to as canonical GRBs in 
DIO) define an subsample t/0095 for the LT correlation. Actually, one has also to consider 
the error bars on the fitted parameters although we remember that b is actually correlated 
to a and being analitically set by Eq. OlOp . When the large error bars are taken into 
account, a can indeed be considered independent on Uth, while the trend with Uth of the 
zeropoint b remains meaningful. We place here a general explanation on the size of the error 
bars presented in Fig. [2l Fig. |5l Fig. [7] and Fig. |8l namely how the uncertainties on the 
calibration parameters have been derived. The size of the uncertainties does not reflect only 
the scatter in the data in the plots, in fact we can note that they are greater than la, since 
they reflect also the intrinsic scatter in the law log L*^ = a log T* + b. Furthermore, in Fig. [S] 
the errorbars represented are not directly obtained from the Equ. |6] but they are computed 
as the median absolute deviation from the median of the parameters. The median values 
of the observables in the GRBs lightcurves present highly scatter since they reflect intrinsic 
i nhomogeneities in the parameters values . As discussed in detail in our previous papers 
(jPainotti et al. 2008t iDainotti et al. 20101 ). the constraints on (a, b, aint) have been obtained 
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Fig. 4.— Same as Fig.[3]but for the UP sample (8 GRBs with u < 0.095). 



by running a Monte Carlo Markov Chain algorithm to explore the parameter space (a, o"mt), 
while b is analytically derived through Eq. lfTOj). To this end, we run two chains, check 
convergence through the Gellman - Rubin test (IGelman fc Rubin 19921 ) and finally merge 
them to estimate the median value and the 68% and 95% confidence ranges. Figs. [3] and 
m show these histograms for the fits to the 77 GRBs with u < A and the U0095 sample. 
Therefore, the error bars in Fig.[2] are determined not only by the scatter in the data, but 
also by the degeneracies among the model parameters. Moreover, being the distributions 
mildly asymmetric, the 68% confidence range should not be meant as la error although 
we use these abuse of terminology for sake of simplicity. When comparing the results of 



^ A caveat is related to the high peak in the b histogram. Because of parameters degeneracy, there wih be 
different possibihties to get a value of b close to the best fit one. Since the code preferentially selects models 
with (a, aint) as close as possible to the best fit parameters, we will get a lot of couples (a, aint) giving almost 
the same b value. In order to show the full b distribution, we have chosen a range much larger than what 
is actually needed so that the central bin (i.e., the one which the best fit b lies within) will be much more 
populated than the other ones hence explaining the peak in the figures. Note also that the first bin in the 
P(— 21ni) plots is less populated because it is the one corresponding to the best fit parameters. In order to 
reach convergence, the MCMC code must first find the best fit and then move away from here so that the 
first bin is not the most populated one. 
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the fits to the different u selected samples, one should therefore compare the histograms on 
(a, 6, aint) and then consider the calibration parameters of different fits in agreement if the 
corresponding histograms well overlap. This is, for instance, the case for the fiducial and 
U0095 samples. We note that the median values of the distributions are different, being 
{a,b,aint) = (—1.04,51.30,0.76) for the fiducial sample and {a,b,aint) = (—1.05,51.40,0.40) 
for the U0095 one. However, the histograms for a well overlap so that we find no statistical 
meaningful difference, while this is the case for both b (although weak) and amt- The 
error bars plotted in Fig. [2] allow a quick check for the samples with varying u making us 
confident that the trends commented above based only on the median values are statistically 
meaningful. 

Up to now, we have interpreted the selection on u as a way to find out the GRBs 
most closely following the WOT model. It is worth wondering whether such a selection 
biases in some way the sample by selecting only, e.g., high luminosity GRBs or the shortest 
ones. To this end, we show in Fig. [5] the median values (with the median deviation) of 
(logL^f , logT^, /3a) as function of the threshold value used for the error parameter u. As 
these plots clearly show, there is no trend of the median values of (log L^, log T*, /?„) with u, 
namely the invariance of u with respect to (log L^^^, log T^) corresponds directly to invariance 
in the samples that result when GRBs are chosen for certain u. Indeed, even neglecting the 
large error bars [^, the median values keep constant showing that the samples selected by 
imposing u < Uth sample the same region in the parameter space (log L^, logT*, /J^). This 
is due to the fact that the definition of u depends directly on cr|. and cr|., . Therefore, the 
direct dependence on the possible biases on the sample depends on the parameter values 
that characterize u. Nevertheless, to be confident that correlation among the parameters 
will not affect the sample selection we have tested also the median values of Pa vs m, because 
there is an indication of a corr elation among logT^ vs especially for the limiting u = 0.095 



sample, (jPainotti et al. 20101 ). All the tests described clearly shows that the selection on u 
is only a way to find out the GRBs following as close as possible the Willingale's model, but 
this criterion don't bias the samples. As a consequence, we can conclude that the smaller 
scatter of the LT correlation for canonical GRBs is not a product of selection effects, but 
rather the outcome of a (still to be understood) physical mechanism. 

A careful inspection of the log L*j^ vs log T* plot suggests that the most deviating points 
are the low luminosity GRBs. We have therefore repeated the above analysis by selecting 



^Note that we are here using the median deviation to characterize the width of the distribution so that, 
strictly speaking, these are not la errors. Note also that the first point in every distribution has a larger 
error bar since the corresponding sample is made out of only 4 GRBs with u < uth and uth < 0.095. We 
still plot this point for completeness although it is likely statistically meaningless. 
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Fig. 5. — Median and median deviation values of (log L^, log T^*, /3a) as function of the 
threshold values of u (left panels) and log (right panel) . 
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Table 1: Results of the calibration procedure for GRBs divided in three equally populated 
redshift bins with {z^in, z^ax) = (0.08,1.56), (1.71,3.08), (3.21,8.26) for bins Zl, Z2, Z3. 



PLT i.^bf: ^bf 1 ij^intjbf^ O'median ^median i^^^int) median 



Zl 


-0.69 


(-1.20, 51.04, 0.98) 


-i-08lS:L^ 


51.05lJ:i 


1 01+0-20 

-'^•'-'-'--0.16 


Z2 


-0.83 


(-0.90, 50.82, 0.43) 


-0.86l°:i« 


50.901°:?^ 




Z3 


-0.63 


(-0.61, 50.14, 0.26) 


-0.581°:}^ 


50.15loi9 


26+°-°^ 



samples with logL^^^ > {log L\)th with {log L*x)th running from 43.50 to 46.50 in steps of 
0.15. As shown in the right panels of Fig.|5l such a selection criterion do not bias the sample 
in (logL^c, logTj^, /3a) hence suggesting that the intrinsic scatter of the LT correlation could 
be reduced by using only moderately bright GBRs. However, since we do not have a physical 
motivation for applying such a criterion, we have not carried out this analysis. 

4. A redshift dependent calibration ? 

The redshift range covered by our GRBs sample is quite large with {zmin^Zmax) = 
(0.08, 8.26) although the distribution is actually quite inhomogenous. Indeed, we have a 
single GRB at z = 8.26 with the second farthest GRBs being aX z = 5.3. Similarly, the 
closest GRBs is at z = 0.08, but the second one is at 2 = 0.12. The redshift distribution 
has played, up to now, no role in our analysis since the LT correlation has been fitted to the 
full GRBs sample thus implicitly assuming that the calibration coefficients (a, 6, (Tint) are the 
same over this wide redshift range. It is worth wondering whether this is actually the case. 
To this end, we have therefore recalibrated the LT correlation dividing the GRBs in three 
equally populated redshift bins. Note that we use here the fiducial sample {u < 4) in order 
to have a good statistics. If we had chosen, for instance, a set with u < 0.3, only 11 GRBs 
per bin would be present thus leading to large errors preventing any comparison. 

The results summarized in Tableland shown in Fig. [6] allow us to draw some interesting 
remarks. Firstly, we note that, although the intrinsic scatter is quite large (mainly because 
of the use of the fiducial rather than the UP sample), the correlation coefficient pir is quite 
large in all the redshift bins thus arguing in favour of the existence of LT correlation at any 
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z. The slopes a for bins Zl and Z2 are consistent within the 68% CL, while this is not 
for bin Z3 where the agreement is present only at the 95% CL level. On the contrary, the 
zeropoint h is consistent among the three bins. Although a trend in the median values of 
a is present, we therefore can not conclude that the LT correlation becomes shallower for 
higher redshift GRBs because of the paucity of the sample and the inclusion of large u GRBs. 
Larger samples with low u values and a more homogenous redshift sampling are needed to 
solve this issue. 

The study of redshift evolution of the LT correlation is particularly interesting in view of 
its application to cosmology. If the calibration parameters had significantly changed with z, 
one could have not used the same set of parameters for all the GRBs as, on the contrary, we 
have usually done. To better clarify this issue, we first remember that the distance modulus 

/x = 25 + 51ogDL(z) (14) 

depends on the cosmological parameters as shown by Eq.(I7]). On the other hand, because of 
the Eq.(|6]), the value of /i for a GRB at redshift z may be inferred by the LT correlation as : 



with the error estimated as : 



25 + -log 



47^fa{Ta,Tp,FaTa){l + zY-'+M 



25 + - <; a log 



T 

-J- a 



l + z 



+ b 



- log [irrUn, Tp, F,T,)(1 + z^-'^^^^ 



(15) 



Dl[z) mlO 

We stress here that the total uncertainty is obtained by adding up the statistical error from 
the propagation of the errors on the involved quantities and the intrinsic scatter. 

We present two tests performed with the aim of understanding if the LT evolves with 
redshift. We denote with fizi and fifid the values of fi estimated from Eq. (fT5|) using (a, b, aint) 
obtained by fitting the Zi and the fiducial samples respectively. Fig. [7] presents the plots 
= fizi — fJ'fid z for the three redshift bins showing that A/x is has a roughly constant 
behaviour in each redshift bin. Moreover, as shown in Fig.|8] (giving instead fizi/f^fid as 
function of z), we can see a flat behaviour too, or at maximum a difference of order ~ 2%, 
much smaller than the typical error bars and hence fully negligible. We therefore argue a 
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1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 

log[rj(l+z)] 



Fig. 6. — \ogL*x,\ogT* correlation divided in the three redshift bins Zl = (0.08,1.56), 
Z2 = (1.71,3.08) and Z3 = (3.21,8.26). With the blue points we have represented Zl 
sample, with the green ones the Z2 sample and with the red points the Z3 sample. The 
respective fitted lines are in the same colours. 
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Fig. 7. — = fizi — fJ'fid as function of the redshift z for the three bins Zl (0.08 < z < 1.56), 
Z2 (1.71 <z< 3.08), Z3 (3.21 <z< 8.26). 
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Fig. 8. — Same as Fig. [7] but for the distance modulus ratio Hzil l^fid- 



(still to be confirmed) redshift dependence of the LT correlation does not preclude its usage 
as a way to construct a GRBs Hubble diagram (ICardone et al. 20091 : ICardone et al. 20091 ) 
for cosmological applications. 



5. Looking for a redshift estimator 

The above analysis has shown that the LT correlation is empirically well motivated and 
not affected by selection effects due to u selection or to redshift dependent calibration. It is 
therefore worth investigating its possible applications as redshift estimator. To this aim, let 
us go back to Eq. (E]) and rearrange it in a different way as follows : 



log[L:^(T,)] = log(47rFx) + 21ogDL(^)-(l-/3a)log(l + z) 

= log {At: Fx) + (1 + /?«) log {l + z) + 2 \ogr{z) + 2 log {c/H^) 



a log 



Tg 

l + z 



(17) 
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where we have denoted with r{z) the integral in Eq.([7]) and, in the last row, we have used 
the LT correlation with the definition of T*. Solving with respect to z, we get : 



(1 + /3, + a) log {l + z)+2 log r(z) = a log + 6 - log {irrFx) - 2 log (c/Ho) (18) 



For the considered cosmological model, the right hand side of Eq. (|T8i) depends only on 
measurable quantities so that one can try solving this equation with respect to z to get an 
estimate of the GRB redshift. There are, however, some preliminary issues that must be con- 
sidered. First, both the observable quantities (T^, Fx, Pa) and the LT calibration parameters 
(a, b) are affected by their own uncertainties. Propagating these errors on the final estimate 
of z is not analytically possible. Moreover, the uncertainties on (a, b) are not symmetric and 
the intrinsic scatter aint (also known with its own asymmetric confidence range) adds to the 
total uncertainty in a nonlinear way. If we denote by Z{p) the solution of Eq.( !T8|) for a 
given set of parameters p = {logT^, log Fx, Pa, o, b}) and neglect the correlations among the 
errors, one should estimate the error on z^st as : 



where the sum runs over the number of parameters. Actually, such a formula can not be 
used since, firstly, we do not have an analytical expression for Z{p) and, secondly, there is 
actually a non negligible correlation among the parameters (for instance, b is determined 
from the value of a and aint)- To fully take into account this issue, for each GRB, we first 
estimate z setting all the observable quantities (logT^, log Fx, Pa) to their central values and 
the calibrationparameter (a, 6) to their best fit values and solve Eq. (fT8|) to get what we 
denote as Zest- □ 

We have applied the above test0 to the fiducial {u < 4.0) and the U0095 samples 
[u < 0.095) finding out that the LT correlation can still not be used as a redshift estimator. 



■^Regarding the uncertainty we can only provide a rough estimate repeating the procedure described for 
a large set of randomly generated values of (a, &), derived by interpolating the (normalized) histograms 
outputted from the Markov chains. We then take the histogram of the z^st values thus obtained to find 
out the 68% confidence range {zrmn,Zmax) finally defining cr{zest) = [{zmax - Zest) + {zest - Zmm)]/'^, i-e. 
we symmetrize the confidence range. It is worth stressing that such an approach implicitly allows us to 
propagate also the error on aint since the zeropoint h is determined from the values of {a, aint) on a case- 
by-case basis. We don't present in the Fig. |9]the errorbars not to cutter the picture and because for the 
reasons discussed above they can be only an approximated estimate of the real error measurements. 

^Remember that, when performing the test, we use the merged chains relative to the fit to the considered 
sample so that, for instance, the best fit (a, b) values are different in the two cases. 
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2 4 6 8 

^obs 

Fig. 9. — Observed vs estimated redshift for the 62 GRBs of the fiducial sample divided in 
three u bin, i.e, u < 0.095 (red points), 0.095 < u < 0.3 (blue points), 0.3 < u < 4 (black 
points) . 

as clearly shown in Fig.O Indeed, defining Az = Zobs — Zest, we get only ~ 20% (28%) 
of GRBs in the fiducial (U0095) sample has \Az/a{zest)\ < 1- Even if we allow for a very 
poor precision considering as acceptable estimates those with \ Az/a{zest)\ < 3, the fraction 
of successful solutions raises only to a modest ~ 53% (~ 57%) for the fiducial (U0095) 
samples. The qualitative agreement about the bad performance of this estimator for both 
the fiducial and U0095 samples is a first evidence that the value of u has no impact on 
the quality of the redshift estimate. This is also shown in Fig.[9] where the points closer 
to the Zobs = Zest line are not the ones with the smaller u values. Actually, such a result 
can be easily understood noting that, as yet demonstrated above, the u selection does not 
bias the samples so that the underlying motivation why this redshift estimator fail applies 
equally to all GRBs notwithstanding how precise is the measurement of (logT*, logL^). 
Actually, the main motivation of the failure is the intrinsic scatter of the points around the 
best fit LT correlation. It is quite easy to qualitatively understand this point considering a 
hypothetical GRB of the U0095 sample with given value of logT^. Because of the intrinsic 
scatter of the LT correlation, its estimated intrinsic luminosity logLx can be off from the 
true one up to aint- Let us denote with L^"^ and L;^* its true and fitted Lx values and let 
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Lx{z) = 'iTTDL{zyFx{l + z)~^^~'^°-\ Our estimate of z is obtained by solving Lx{zest) = L^x^, 
while it is Lx{zobs) = L^x^^ ■ Should L^x be larger than L^"^, we should increase Dl{z) to 
compensate for the difference thus overestimating the redshift {zest > Zobs) with the opposite 
effect for the L^^ < ^x"' 

case. This is indeed what happens in our case. Looking at 
the residuals of the LT correlation, we find that Lx is, on average, underestimated (i.e., 
L^^ < L^"*^) for the higher redshift GRBs so that we expect Zest < Zobs which is indeed what 
we find looking at the points with z > 2.5 in Fig.[9l while the opposite result takes place for 
very low z GRBs. We can therefore conclude that the LT correlation works satisfactorily 
well as redshift estimator only for GRBs in the redshift range (0.5, 2.0), while gives severely 
biased results for smaller and larger z GRBs. 

Although these results are quite discouraging, it is nevertheless worth wondering whether 
the situation can be improved with future data. To this end, we have simulated a U0095 sam- 
ple generating (logT*, I3a, -z) values from a parent distribution closely mimicking the observed 
one for the fiducial samplqf|. We then set logL^ extracting from a Gaussian distribution 
centered on the value predicted by the LT correlation and with width equal to the intrinsic 
scatter. We then use these values to estimate Fx and add noise to all quantities so that the 
relative errors are the same order as the present day ones. We generate A/" GRBs and fit 
them with the same procedure adopted to find the LT calibration coefficients and use these 
fake Markov chains as input to the redshift estimate procedure. 

It turns out that increasing the sample is not a useful way to improve the performance of 
the redshift estimator. Indeed, we have found that, with J\f ~ 50, the fraction of GRBs with 
\Az/a{zest)\ < 1 first increases to ~ 34% and then decreases to ~ 20% for M ~ 200, while 
(Az/zobs) — —17% for both ~ 50 and J\f ~ 200, a significant improvement with respect to 
the value quoted above, but still not fully satisfactory considering that rms{Az/zobs) — 75%. 
It is somewhat counterintuitive that increasing the sample does not improve the quality of the 
redshift estimator. Actually, such a result could be anticipated noting that a larger sample 
leads to stronger constraints on the (a, b, (Tint) values, but do not change the intrinsic scatter 
which is the main source of possible mismatches between the true and fitted GRB luminosity. 
Motivated by this consideration, we therefore perform a second test artificially lowering the 
intrinsic scatter aint, but setting the best fit (a, 6) parameters to those derived from the 
fit to the U0095 sample. Indeed, for A^ = 50 and aint = 0.20, we get {Az/zobs) — —3% 
with rms{Az / Zobs) ^ 28% and 46% (87%) of GRBs with \Az/a{zest) < 1(< 3). Again 
increasing the sample to A/" ^ 200 has not a significant impact, while a stronger impact 



^Note that this choice is motivated by the poor statistics of the present U0095 sample. We have, however, 
checked that the U0095 GRBs cover the same region in the (log Ta, /3a, -z) parameter space as the fiducial 
ones. 
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is obtained setting amt = 0.10 giving {Az/zobs) — —0.6% with rms{Az / Zobs) — 16% and 
f[\Az/zobs\ < 1) — 66%. These resuhs convincingly show that the LT correlation could be 
used as a redshift estimator only if a subsample of the canonical GRBs could be identified 
in such a way to reduce the intrinsic scatter to aint = 0.10 — 0.20. It is worth wondering 
whether assembling such a sample is indeed possible. Actually, a detailed answer can not be 
given since our U0095 sample is too small to find out some indicator that can help to find out 
GRBs less scattering from the best fit LT correlation. A visual inspection of the fit residuals 
makes us roughly argue that a^nt could be reduced using only 5 out of the 8 U0095 GRBs 
which represent ~ 8% of the fiducial GRBs sample. If we assume this fraction as a constant, 
one should assemble a sample of ~ 600 GRBs with measured values of (logT*, log L^, /3a, z) 
to get ~ 50 GRBs to calibrate the LT correlation with aint ~ 0.20. While this is for sure an 
ambitious task, it is worth noting that it is still possible that a smaller sample is enough to 
find out the observable properties of these GRBs thus allowing an easier search and reducing 
the number of GRBs to be followed up for the z estimate. 



6. Summary 

The analysis presented here have shown that the LT correlation, for both the fiducial 
and UE samples, is not affected by selection effects induced by the u threshold selection or by 
the implicit assumption of redshift independence of the calibration parameters. In particular, 
the selection on u does not bias the distribution of the (logL^^ , logT^, /3a) quantities thus 
showing that the canonical GRBs {u < 0.095) in the t/0095 sample are indeed distributed 
preferentially in the upper part of the LT plane. This is a further evidence that the afterglow 
light curves which are smooth and well fitted by the W07 model indeed define a physically 
homogenous class with the remarkable feature of obeying a well defined empirical correlation. 
Furthermore, the analysis presented also pinpoints the existence of a well defined correlation 
of the X-ray spectral index Pa with the rest frame break time T* which deserves further 
analysis. 

As an important result, we have also shown that, although a shallowing of the LT 
correlation for higher z GRBs can still not be totally excluded, its impact on the distance 
modulus estimate is neg hgible thus vahdating the usage of this correlation as a new indepen- 



dent cosmological tool (jCardone et al. 20091 ). As a first application, Cardone et al. (2010) 
have indeed derived the Hubble diagram using the LT correlation only and shown that, when 
combined with other distance probes (such as Type la Supernovae and Baryon Acoustic Os- 
cillations), GRBs are a valuable tool to constrain the cosmological parameters. To this end, 
Cardone et al. (2010) have used the full fiducial sample to increase the statistics, but at the 
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price of including GRBs with large errors on the distance modulus. It is worth wondering 
how large a GRB sample should be to improve the constrain ts on the cosmologic al parame- 
ters. Such a problem has yet been addressed by some of us (jCardone et al. 2009 ) using the 
Fisher matrix analysis and the first version of the LT correlation. There, we have shown 
that combining a sample of 200 GRBs with a SNAP - like SNela sample allows to determine 
the matter density parameter Qm and the dark energy equation of state parameters (wq, Wa) 
within 0.019, 0.036, 0.020, respectively. In particular, GRBs are of extremely importance to 
constrain Qm giving an improvement in precision of a factor 4 with respect to the case when 
SNela only are used. Although these results refer to the first version of the LT correlation 
and thus refer to a sample with no selection on u, they qualitatively hold also in our case 
since the basic inputs to the Fisher matrix analysis are essentially the same. Actually, having 
made no cut on u, the quoted results are likely to be quite conservative since the u selection 
allows to reduce the scatter and hence the error on the distance modulus thus increasing the 
efficiency of GRBs with respect to the case considered in Cardone et al. (2009). 

We have, finally, investigated the possibility to use the LT correlation as a redshift 
estimator obtaining actually discouraging results for both the fiducial and U0095 samples. 
Having qualitatively discovered the reason of this failure, we have shown that reducing the 
intrinsic scatter of the LT correlation could help to calibrate an improved LT correlation that 
could work as a tool to estimate the GRB redshift from the analysis of its X - ray afterglow 
lightcurve. However, we are not sure if it is possible to obtain a reduced intrisic scatter of 
the correlation with the real data measurements related to the U0095 sample. 
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